Anomalous inter-layer atomic transport and the low-temperature amplification of 

surface instability in Al(lll) 
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Anomalous inter-layer atomic transport of deposited impurity atoms on Al(lll) has been found by 
constant temperature molecular dynamics simulations. The low-energy deposition of Pt on Al(lll) 
leads to oscillatory adsorbate-substrate interaction and to a low-temperature ballistic injection of 
the deposited particles to below the topmost layer. The ultrafast injection of a Pt atom coincides 
with the ejection of a substrate atom to the surface (ballistic replacement or exchange mechanism) . 
This is in agreement with the experimental findings in which thin film rich in Al has been found on 
the surface after deposition of few MLs of Pt. We attribute the anomalous inter- layer transport to 
the size- mismatched impurity/host interaction and we point out the role of atomic size mismatch 
in biasing towards intermixing. The deposition induced low-temperature disordering of surface Al 
atoms with few transient Al adatoms has also been observed. The ultrafast injection of the impurity 
particles to the substrate is assisted by the transient out-of-plane and lateral circulating motion of 
few surface atoms arranged nearly in a hexagonal symmetry. The atomic injection of Pt can be 
regarded as a superexchange process driven by the transient and collective motion of few surface 
atoms. A chaotic surface state assisted mechanism could also be a possible explanation of the 
unexpected phenomenon. 
PACS numbers: 68.43. Jk, 68.35.Fx, 66.30.-h 
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I. INTRODUCTION 

The understanding of the fundamental atomistic trans- 
port processes leading to the formation and morphologi- 
cal evolution of nanoscale surface features is lacking [1]. 
Many classical macroscopic (continuum and mesoscopic) 
models for diffusion and morphological evolution lose 
their validity in the nanoscale [1,2]. Therefore the atom- 
istic level understanding of the driving forces and mecha- 
nisms governing atomic transport involved in the synthe- 
sis and organisation of nanoscale features in solid state 
materials including atomic intermixing, is inevitable [1]. 

The most well known atomic diffusion mechanisms 
of adatoms on solid surfaces are the hopping diffusion 
on the surface from one site to another over a bridge 
site and the atomic site exchange when an adatom en- 
ters the topmost layer and a surface atom becomes an 
adatom [3]. The atoms deposited at surfaces could un- 
dergo, however, few exotic (or anomalous) transport pro- 
cesses such as random walk or Levy flight [3-5] , Friedel- 
type adsorbate-adsorbate oscillatory interaction driven 
relaxations [3,6-9], intermixing and surface alloying be- 
tween immiscible elements [7,10,11,14], long range sur- 
face mediated lateral mass transport within the topmost 
layer [12] and very fast kinetics with low diffusion bar- 
rier driven by quantum size effects [13]. The low-energy 
ion-bombardment induced ejection of Al atoms to the 
Al(lll) surface results in adatoms whith ultrafast mo- 
tion on the surface leading to the formation of nanodots 
[15]. These atomic migrations can be characterized by 
the lateral transport of atoms on the surface or whitin 



the topmost layer. Inter-layer transport occurs only be- 
tween the surface and the topmost layer or at step edges, 
however, the most of these processes can not usually be 
considered as anomalous, e.g. the rate of the process can 
be described within the framework of the Arrhenius equa- 
tion. Anomalous atomic processes are athermal and the 
mean free path of the atomic jumps can reach the few 
times of the nearest neighbor distance [3,4]. Also, the 
mean square of atomic displacements scales nonlinearily 
in the course of time of diffusion [4,16] or the nanoscale 
growth rate of interfaces follows linear growth kinetics 
(nonparabolic growth) [2] . 

However, only few examples are known for inter-layer 
mass transport (atomic jumps normal to the surface) 
which can be considered as anomalous. Even in these 
articles the authors not always realized the anomalous 
nature of intcr-layer transport. The assistance of instan- 
tenously high kinetic energy to inter-layer transport dur- 
ing the deposition of Au on Ag surface has been reported 
[17]. Interdiffusive Stranski-Krastanov growth mode in 
an analysis of scanning tunneling microscopy data and 
the burrowing of Au particles in Ag(llO) when anneal- 
ing has been applied during simulations (replacement 
diffusion mechanism) have also been found [18]. In an- 
other study a complicated exchange mechanism is shown 
with ballisticaly injected deposited Au on Ag(lll) [19]. 
Sprague et al. has been published interface mixing on 
impact (athermal mixing) for the case of Pt deposition 
on Cu and which is explained by the attractive potentials 
of the substrate atoms and a large thermodynamic bias 
toward film-substrate mixing [20]. The mechanism of 
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the seemingly different transient atomic migrations could 
hopefully be explained in a common framework. 

The understanding of how adsorbed impurities and 
growing film constituents affect interdiffusion to the sub- 
strate needs further studies of the atomistic mechanism 
of intermixing. For instance, using vapor deposition of 
various transition metals on Al, it has been found that 
the intermixing (IM) length is anomalously large in cer- 
tain cases [21]. It has also been characterized that IM is 
not driven by bulk diffusion parameters [21] nor by bulk 
thermochemistry (such as heats of alloying) during ion- 
beam mixing [22,23]. The segregation of Al during vapor 
deposition of Pt on polycrystalline Al together with the 
fast reactive diffusion of Pt to the bulk have also been 
found [21,24]. In another study the adsorption of Pt on 
polycrystalline Al leads to the formation of surface alloys 
[25] which are rich in Al [26]. 

Anharmonic effects, surface disordering and premelt- 
ing has also been studied intensively at fee (111) metal 
surfaces [27,28]. It has also been concluded that surface 
thermal expansion and lattice dynamics might be anoma- 
lous of the (111) surface layer of Al [27,36] hence Al can 
be described by unusual surface instability. It would be 
interesting if e.g. the surface of Al would be subjected 
to small external perturbation such as the deposition of 
an impurity atom. Owing to the anharmonic surface be- 
havior of Al(lll) the amplification of surface instability 
could be expected. 

A number of unconventional atomic transport phenom- 
ena have been explained by atomic size- mismatch (ASM) 
of the constituents [10,16,29 35]. Surface alloying and 
intermixing between many immiscible and even misci- 
ble elements can not solely be explained by thermody- 
namic and chemical forces such as heat of alloying and 
mixing or cohesive energies [10,11,14,22,23]. The puz- 
zling nature of intermixing (IM) might be due to still an 
unknown and an unconventional mechanism which has 
not been explored yet or overlooked in the last decades. 
In our previous reports we explained IM as an interfa- 
cial anisotropy driven process in diffusion couples and 
pointed out the role of mass and size-anisotropy during 
ion-bombardment induced interdiffusion [15,16,22]. 

In the present article, we show that anomalous 
atomic transport occurs during vapor deposition in size- 
mismatched systems, such as during the thin film growth 
of Pt and other transition metals on Al(lll). The inter- 
diffusion of the deposited particle is driven by an oscilla- 
tory interaction with the surface atoms of the substrate. 
Moreover, we point out that the anomalous character of 
the atomic transport of the deposit can be tuned by ad- 
justing the ASM of the constituents. We present an ex- 
change mechanism with a concerted motion of few surface 
atoms for entering the uppermost layer by the deposited 
particle. Also, we show that large atomic vibrational dis- 
placements and anharmonic effects might occur not only 
at high-temperatures close to the melting point, but also 
at ultra-low temperatures induced by impurities at the 
surface. 



II. THE SIMULATION METHOD 

A. General properties 

We give the detailed outline of the employed simula- 
tion approach. Classical constant volume tight-binding 
molecular dynamics simulations [43] were used to sim- 
ulate soft landing and vapor deposition of Pt atoms on 
Al(lll) and also on other substrates (such as Cu(lll)) 
surface at - K using the PARC AS code [37]. The 
PARCAS MD code has widely been used for the study of 
various atomic transport phenomena in the last few years 
[22,38]. A variable timestep and the Berendsen temper- 
ature control is used [39,40]. The simulation uses the 
Gear's predictor-corrector algorithm to calculate atomic 
trajectories [39]. The maximum time step of 0.05 fs is 
used. We consider the coupling of our simulation cell to 
a heat bath by inserting stochastic and friction terms to 
the equation of motion yielding a Langevin type equation 
(see details in refs. [39-41]). The equation of motion can 
be written as follows, 

rriiVi = Fi - m.i'jiVi + R{t), (1) 

where Fi is the force on atom i and R{t) is a Gaussian 
stochastic variable. The damping constants 7j determine 
the strength of coupling to the thermostat. The global 
coupling to the heat bath can be adjusted by the so called 
Berendsen temperature which we set to 70 K. The system 
couples to a heat bath not only globally via the damp- 
ing constant but is also locally subjected to random noise. 
Stochastic forces and damping applied at the cell borders 
(lateral x — y boundaries as heat sink) of the simulation 
cell to maintain constant temperature conditions and the 
thermal equilibrium of the entire system (coupling to the 
heat bath) . The top of the simulation cell is left free (the 
free surface) for the deposition of Pt atoms. The bot- 
tom layers are held fixed in order to avoid the rotation 
of the cell, system. For simulating deposition it is ap- 
propriate to use temperature control at the cell borders. 
This is because it is physically correct that potential en- 
ergy becomes kinetic energy on impact, i.e. heats the 
lattice. This heating should be allowed to dissipate nat- 
urally, which means temperature control should not be 
used at the impact point. Periodic boundary conditions 
are imposed laterarily. The observed anomalous trans- 
port processes are also observed without periodic bound- 
ary conditions and temperature control. Further details 
are given in [37] and details specific to the current system 
in recent communications [22,23]. 

The size of the simulation cell is 80 x 80 x 42 in- 
cluding 16128 atoms (with a fee lattice). 15 active MLs 
are supported on 3 fixed bottom monolayers (MLs). We 
find no dependence of the anomalous atomic transport 
properties of the deposited atoms on the finite size of the 
simulation cell. Deposited atoms were initialized normal 
to the (111) surface with randomly selected lateral posi- 
tions 4 — 5 A above the surface. 
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The kinetic energy of the deposited particles are ^0.1 
eV or in the case of uhrasoft landing it is nearly zero 
eV. The impurity (deposited) particle is placed 4 — 5 A 
above the (111) surface of the substrate. No acceleration 
of the deposited particle is observed above 6 A the sur- 
face. We also analyze the acceleration of the deposited 
particles and calculate the arrival energy at the surface 
of the substrate. Using the history (movie) file we col- 
lect the atomic positions of the moving substrate atoms 
which have kinetic energy above a certain threshold value 
in order to get the pattern of atomic trajectories during 
the deposition events. At K this value is ^ 0.01 eV. In 
order to make a statistics of impact events we generated 
100 events with randomly varied impact positions. 

In order to follow the time evolution of atomic motions 
the mean square of the atomic displacements (R^) (MSD) 

has been calculated. (R^) — J^fi^ii^) ~ = 0)]^, as 
obtained by molecular dynamics simulations and (ri(t) 
is the position vector of atom 'i' at time t) and N is 
the number of atoms included in the sum, respectively. 
We also calculate (i?^) for the impurity atom only, hence 
in this case (i?^) (ri(i) — ri{t = 0))^. In other cases 
lateral or vertical (out-of-plane) components are included 
in (i?^) for the substrate atoms. We folow the evolution 
of (J?^) during few events in order to get a reasonable 
statistics of atomic motions. 



B. The interaction potentials 

We use the many-body Cleri-Rosato (CR) tight- 
binding second-moment approximation (TB-SMA) inter- 
action potential to describe interatomic interactions [43] . 
The CR potential is formally analogous to the embedded 
atomic method (EAM, [52]) formalism, e.g. the poten- 
tial energy of an atom is given as a sum of repulsive pair 
potentials for the neighboring atoms (usually for the first 
or second neighbors and a cutoff is imposed out of this 
region) and an embedding energy that is a function of 
the local electron density given as follows [52] , 



(2) 



where is the distance between atoms i and j and its 
neighbors. There are many functional forms are available 
for the density pi and for the embedding function F[pi] 
[52]. We utilize EAM functional forms in the code for 
F[pi] and for the density p similar to that given in refs. 
[52,53]. The EAM routine in the code employs a cubic 
spline interpolation for the evaluation of the EAM poten- 
tials and their derivatives (forces) starting from various 
kind of input potentials given in discrete points as a func- 
tion of rij (the number of points per functions is 5000 in 
this study). 

Within the second moment tight-binding approach, the 
band energy (the attractive part of the potential) reads. 
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where Tc is the cutoff radius of the interaction and ro is 
the first neighbor distance (atomic size parameter). 

The repulsive term is a Born-Mayer type phenomeno- 
logical core-repulsion term: 
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The parameters {^,q,A,p,ro) are fitted to experimental 

values of the cohesive energy, the lattice parameter, the 
bulk modulus and the elastic constants Cn, Ci2 and C44 
[43] . The summation over j is extended up to fifth neigh- 
bors for fee structures [43]. The total cohesive energy of 
the system is 
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where U{r) is the total potential energy of the entire sys- 
tem as a function of the space coordinate (radius vector) 
r since the system is energy conservative, the space inte- 
gral over the Newtonian interatomic forces gives the total 
energy of the simulation cell, r can be replaced by the 
internuclear separation rij in the pair interaction term 
V{rij) and by the position vector of the electron den- 
sity in the density function Pi{r) in Eq. (2). Formally 
we give all the terms in Eqs. (3)-(4) in the same way 
as for EAM potentials except that F[pi] is calculated in 
a tight-binding form as given in Eq. (3). We prefer to 
use the CR potential because tq can be tuned in this 
case which is proportional to the atomic size-mismatch 
in the heteronuclear potential (the proportionality will 
be discussed later in section III.). Using the CR poten- 
tial we consider the interaction between two atoms and 
the interaction with their local environment usually up 
to the second neighbors. Out of this region a potential 
cutoff is imposed. The cutoff radius Vc is taken as the sec- 
ond neighbor distance for all the interactions. We tested 
the Al-Al and the Al-Pt potential at cutoff radius with 
third and larger neighbor distances and found no consid- 
erable change in the results. The CR potential gives the 
physically reasonable representation of metals and com- 
putationally is efficient hence nanoscale atomic clusters 
with large number of atoms (up to millions of atoms) can 
be simulated. This type of a potential gives a very good 
description of lattice vacancies, including migration prop- 
erties and a reasonable description of solid surfaces and 
melting [43] . Since the present work is mostly associated 
with the elastic properties, melting behaviors, surface, in- 
terface and migration energies, we believe the model used 
should be suitable for this study. In a recent report we 
found that the CR potential describes properly Al and 
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Al/Pt [15] and provides reasonable results for the ion- 
bombardment of Al(lll) in agreement with the experi- 
mental results [3]. The CR potential correctly provides 
the adatom binding and dimerization energies, adatom 
island formation upon ion-bombardment in agreement 
with experiment [56]. Recently it has also been shown, 
that the CR potential remarkably well describes diffusion 
in liquid Al [45,46] and energetic deposition of Al clus- 
ters on Al [47]. In order to be more convincing in the 
accuracy of the TB-SMA model we also used another 
parameterization set for Eqs. (3)- (4) obtained by first 
principles augmented-plane-wave calculations [44]. This 
parameterization offers a satisfactory agreement with the 
experimental data for thermal expansion coefBcient, the 
temperature dependence of the atomic mean-square of 
displacement and the phonon density of states of com- 
pounds [44]. However, we find no serious difference in 
the final results obtained for deposition between the CR 
and the parameterization of Papanicolau et al.. There- 
fore we use the original CR parameterization set for Al. 

Despite the empirical nature of the EAM and CR po- 
tentials their accuracy seems to be sufficient for studying 
materials under nonequilibrium conditions such as anhar- 
monic effects and thermal expansion of various fee met- 
als [27,28]. In the present paper we use these potentials 
for the investigation of low-temperature disordering pro- 
cesses with transient surface atomic vibrations for which 
the EAM and CR potentials should be adequate [27,15]. 

For the crosspotential of substrate atoms and Pt we 
employ an interpolation scheme [15,22,48] using the ge- 
ometrical mean of the elemental energy constants and 
the harmonic mean for the screening length of Eqs. (3)- 
(4) are taken as in refs. [48,49]. The CR elemental po- 
tentials and the interpolation scheme for heteronuclear 
interactions have widely been used for MD simulations 
[11,14,22,31,50]. Recently CR interpolated crosspotential 
has also successfully been used for Ti/Pt in agreement 
with our experimental results [16]. The scaling factor ro 
(the heteronuclear first neighbor distance) is calculated 
as the average of the elemental first neighbor distances. 
The AlPt potential is fitted to the measured effective 
heat of mixing in the cubic AlPt (AiJ —100 kJ/mol) 
with a melting temperature of 1870 K [51] which is re- 
produced by our Cleri-Rosato crosspotential within the 
range of 1800 ± 100 K. In order to adjust AiJ in the Al- 
Pt potential (which is proportional to the strength of the 
interaction and to the heat of alloying in the AlPt alloy) 
the preexponential parameter ^ in Eq (3) is set to ^ ~ 3.0 
[15]. Adjusting ^ in Eq (3) one can tune the depth of the 
crosspotential well which is proportional to A7J [15,23]. 

We find, however, that the interdiffusive features of 
Pt in Al does not depend significantly on heat of mixing 
built in the potential in accordance with our earlier find- 
ings for ion-bombardment induced intermixing in Ti/Pt 
[23]. The injective (ballistic) mixing of Pt takes place on 
a broad range of AH values including purely repulsive 
crosspotential (Aif = 0, ^ « 0, attractive term given in 
Eq. (3) is cancelled). This is because anomalous diffusion 



might not be driven by thermodynamic bias (athermal 
mixing) [22]. The real driving force of athermal (bal- 
listic) mixing is far from being understood clearly yet 
[22,23]. 

Depositing Pt on Al(lll) using Ercolesi-Adams (EA) 
Al-Al potential based on the embedded atomic method 
[52,54] we get also injection when the cross pair-potential 
is obtained by the Johnson's scheme [55] which reads as 



VAB{r) = 



^F4r) + ^yB(r) 

PA PB 



(6) 



Hence the heteronuclear pair potential is given as the 
function of the elemental potentials. To avoid singular- 
ities, cutoff must be imposed for the density functions 
Pa and pb equal or greater than the cutoff distance for 
the VA{r) and Vsir) potentials. The cross-embedding 
function is given as the elemental average of F[pi\ func- 
tions. The EA potential has been fitted to ah initio 
atomic forces of many different configurations using the 
force matching method [54] . Recently, it has been shown, 
that the EA potential correctly describes thermal expan- 
sion behavior [54] and anharmonic effects at Al surfaces 
[36,27]. Hence the transferability of the EAM potentials 
to nonequilibrium conditions of Al seems to be suflacient. 



III. THE ATOMIC SIZE-MISMATCH CONCEPT 

Since we study the dependence of atomic dynamics on 

ASM we briefly discuss our approach to ASM. To ac- 
count for ASM the lattice mismatch (LM) concept has 
widely been used in the literature which is given as the 
ratio of the lattice constants of the film and substrate 
constituents [1,10,16,29 35]. In the present paper we use 
another quantity to describe size-mismatch in diffusion 
couples. The motivation for this is that LM seems not to 
explain interdiffusion and surface alloying in a number of 
binary systems [21]. Also, the observed asymmetry of in- 
termixing in various metal film/substrate couples can not 
be understood within the LM concept [16,21]. In the LM 
concept it is assumed that intermixing and many other 
properties of binary systems (e.g. growth modes) are 
primarily determined by the ratio of atomic sizes (lattice 
constants) of the pure elemental phases. Instead we in- 
troduce the quantity atomic size mismatch (ASM) which 
is the ratio of the first neighbor distance tq of the de- 
posited impurity in its alloy phase with the substrate 
(rQ™'^, AlPt alloy in this case) to the ro of the substrate 
in its pure elemental phase (rg). If no alloy exist between 
the constituents than Tq™'" is simply the average of the 
elemental rp values. Hence in the ASM concept the ratio 
5 ASM = ''o™ */''o depends on the cross-interaction be- 
tween the impurity and substrate atoms. It is reasonable 
because the impurity atom when deposited on the sub- 
strate interacts directly with the substrate atoms hence 
mainly the heteronuclear interaction determines ASM at 
low impurity concentration and not its pure elemental 
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homonuclear r™ value. With increasing film thickness 
and film coverage the lattice constant of the film will be 
more and more important ingredient of the overall lattice 
misfit. In other words within this picture of ASM the 
atomic volume (size) is proportional to the local ro value 
which is strongly interaction and environment dependent. 
The same impurity can be described by different ro value 
in its pure and in an intermixed environment. In a strict 
sense the average ro = (rg™ + rQ)/2 is the correct descrip- 
tion of the heteronuclear first neighbor distance in the 1:1 
alloy phases. In nonstochiometric alloys or in inhomoge- 
neous phases (such as intermixed nonequilibrium phases) 
ro™'* fluctuates around the elemental average and locally 
rQ™'* can reach different values depending on the local 
stochiometricity of the system. Nevertheless, the average 
ro seems to be a reasonable approximation for interfacial 
mixing [15,16]. 

In Pt/Al the magnitude of intermixing is not sensitive 
to the choice of ro for the deposition of Pt on Al within 
the range of [2.8; 3.0] A, hence the elemental average of 
r^^'" w 2.85 A is a reasonable choice. We prefer to use in 
most of the presented results the CR potential because 
the first neighbor distance ro can be tuned in the het- 
eronuclear potential in this case which is proportional to 
ASM. Using EAM potentials such as given in Eq. (6) no 
effect of ro can be studied. 

IV. RESULTS AND DISCUSSION 

In the rest of the article we present results which seem 

to be rather surprising and unexpected. Since we find 
no reason on the methodological side to discard them 
(the employed simulation approach and the interaction 

potentials should be adequate for the problem) we feel it 
important to share these atomistic results with a wider 
community. 

The simulation of vapor deposition leads to an unex- 
pected result. The deposited atoms, independently of the 
energy of deposition, intermix with ultrafast atomic ex- 
change entering the top Al(lll) layer even at K within 
a ps if initialized < 4 — 5 A far from the surface. Also, 
the particle not only enter but becomes and interstitial 
atom between the two upper layers. The acceleration of 
the deposited particles have been found starting at ~ 5 
A distance from the surface. The injection of Pt goes 
together with the ejection of one Al atom to the surface. 
The deposition of 1 ML of Pt leads to the formation of 
an adlayer rich in Al in agreement with the experimen- 
tal findings [21,24]. The ultrafast interdiffusion of Pt to 
Al is rather surprising, and similar superdiffusion during 
soft landing of atoms on solid surfaces has been reported 
only for Pt on Cu(lll) [20]. The animation of the atomic 
injection can be seen in a web page [42]. The unexpected 
mechanism of Pt interdiffusion might not be an artifact 
of the CR potential. We have tested recently the various 
surface features of CR Al potential [15] and got fairly 



nice results. Hence, wc do not attribute the anomalous 
behavior of Pt soft landing on Al(lll) to the artifact of 
the CR potential. We find direct injection only in the 
case of certain transition metal elements around Pt in 
the periodic table, such as Ir, Au. Other elements, such 
as Cu intermixes to Al after the deposition of dozens of 
atoms. Moreover, we do believe that the fast reactive 
interdiffusion reported in Pt/Al and in other diffusion 
couples must be due to anomalous diffusion which is over- 
looked in the literature. We get also a rapidly increasing 
intermixing length during the simulation of vapor deposi- 
tion [42] , although still we are far from the experimental 
value of ^ 50 A [21] (results are not shown in this paper). 
We deposit few MLs of Pt which is insufficient to reach 
this value. In this article, however, it is not OTir inten- 
tion to reproduce the experimental intermixing length. 
We rather focus on the understanding of the elementary 
atomic migration step of interdiffusion. We will examine 
in the rest of the paper the details of the mechanism of 
the ballistic interdiffusion of Pt. 

Although we study the mechanism of an atomic jump, 
we discuss the details of an atomic transport step as a 
first step of the reactive diffusion to the bulk reported 
in many papers [21,24]. A number of such single atomic 
jumps leads to a diffusional process through the atomic 
layers of the substrate leading to the large intermix- 
ing length reported for Al [21]. We find the nonlinear 
(paraboHc) scaling of (i?^) as a function of time dur- 
ing interdiffusion which suggests that the atomic inter- 
layer transport is anomalous [4]. An atomic transport 
is anomalous if the exponent a > 1 or a < 1 in the 
expression (i?^) oc t", where, t is the time variable [4]. 

We also vary the first neighbor distance ro in the het- 
eronuclear potential, in order to point out the sensitivity 
of interdiffusion to ASM. The ro of the heteronuclear 
interaction is constructed as the average of the elemen- 
tal values [22,48]. As we have already pointed out, the 
size-mismatch of the binary system is proportional to the 
ratio of 5 asm = r^'^*/r|^'. In this particular case the Pt 
atom is the smaller particle (ro « 2.8) and ro « 2.9 in 
Al. Moreover Al has a highly anharmonic homonuclear 
interaction potential (see upper Fig 1). Due to the an- 
harmonic behavior of Al there is a tendency for increased 
Al-Al distances (r^' > 2.9) on the surface (see upper Fig 
1). If we set ro ~ 3.1 A, no injection is foimd (This is 
because Sasm ~ 1)- Ballistic atomic injection can be 
suppressed by tuning ASM. We also test deposition in 
other size-mismatched systems, such as Ni/Au. Injec- 
tion occurs in this system, when unphysical values are 
used (ro < 2.0 A). 

In order to understand the details of the injective 
anomalous atomic transport the shape of the deposit- 
surface interaction energy potential is calculated. In 
lower Fig 1 we show the binding energy (potential en- 
ergy) of the Pt atom as a function of the distance from 
the surface. An oscillatory impurity-substrate interac- 
tion potential is found. The amplitude of the oscillation 
is large (few eVs). The oscillation of the potential is 
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FIG. 1. Upper panel: The profiles of tlie potential energy 
as a function of the internuclear separation (rjj , A, for Al, 
Pt and for AlPt using the CR = 3.0 in Eq. (3)) and 
Ercolcsi- Adams EAM (EA) potentials ( Al only) . Lower panel: 
The oscillatory potential energy profile of Pt (straight line) as 
a function of the distance («) of the Pt atoms from the surface 
of Al(lll) (the position of the surface is at A, ro ~ 2.83 A). 
The potential (binding) energy of Pt with ro « 3.1 A is given 
with a dotted line. 



caused by the transient out-of-plane movements of sur- 
face Al atoms. We see also the reflection of the cap- 
tured Pt atom from the second Al monolayer (at z « —2 
A). Oscillatory interaction profiles have already been 
reported on the surface of metals between adatoms [9], 
however, no reports have been found for interdiffusion. 

When the Pt approaches the surface one or two Al 
atoms are released from the surface which, however, im- 
mediately return back to the top layer. We show the 
enlarged animation of the active region during injection, 
which can be find at a web page [42]. It is more or less 
clear from this animation that the release of few Al atoms 
to the surface induces the injection of Pt. Finally, the in- 
jection of the Pt goes together with the ejection of an Al 
adatom and the Pt atom is captured below the top layer. 
Prior to the injection we can see a specific mechanism. 
Two or three Al atoms are released ballisticaly towards 
the approaching Pt atom and return back to the sur- 
face in less then 0.1 ps (this Al atom "receives" the host 
atom) . The maximum vertical amplitude of the transient 
Al atoms can reach ~ 2 A. This ultrafast out-of-plane 
vibration of the surface Al atoms do not lead to injection 
yet. Typical This strongly anharmonic behavior of the 
Al surface is shown in the upper panel of Fig. 2 where 
the trajectories of the transient vertical jump of few Al 
atoms to the surface (adlayer) can be seen. The sharp 
potential energy wells in lower panel of Fig 1 also occur 
due to the transient out-of-plane processes. 

Haftel et al. reported an unconventional exchange 
mechanism for Au/Ag(lll) diffusion process involving 
multiple substrate atoms which catalyze interdiffusion 
[19]. They used, however, a reversed annealing simu- 
lation technique which speeds up atomic jumps. Nev- 
erthclcs, they could conclude that the activation energy 
of this specific mechanism is very low. Moreover a col- 
lective double exchange process has been found which 
can be activated ballisticaly by an incoming deposited 
atom. These features arc very similar to that found for 
the Pt/Al system. However, we repeated the simulation 
for Au/Ag(lll) (without reversed annealing) at nearly 
K using the same simulation conditions outlined for 
Pt/Al, and foimd no injection of Au to Ag(lll). Also, 
at elevated temperatures (simulations have been carried 
out up to 300 K) the system does not show up any in- 
tcrdiffusive behaviors up to few picoseconds (not even 
on a longer time scale). However, if we set ro < 2.8 A 
(if we take the elemental average, ro — 2.9 A), Au does 
intermix to the Ag(lll) topmost layer within 2 ps. 

Hence few diffusion couples are only available until 
now, including Pt/Cu [20] and the Pt/Al(lll) (also 
Ir,Au/Al(lll)) systems presented in this work which 
show apparent anomalous intermixing behavior. Few 
other systems, such as the Au/Ag(lll) couple could be 
the subject of a classification of anomalously interdiffus- 
ing systems. Under the effect of forced conditions, such 
as the low-energy ion-sputtering of the Pt/Ti interface 
we also find super-interdiffusion (ballistic jumps) of Pt 
in the Ti substrate and also the nonlinear scaling of {R^) 
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has been obtained from simulations [16]. The anoma- 
lous nature of interdiffusion in Pt /Ti has also been con- 
firmed by Auger electron spectroscopy depth profiling as 
a well-defined long range tail in the concentration profile 
of Pt [16]. This finding provides further evidence for the 
anomalous mass transport behavior of Pt in various me- 
dia. In most of the diffusion couples with no interfacial 
anisotropy (atomic mass and size isotropy) no anoma- 
lous features have been found [16] and vanishingly small 
intermixing length have been measured and calculated 
[22]. However, we do believe that anomalous intermixing 
and impurity caused disordering of the surface of vari- 
ous solids are more widespread in nature than previously 
thought. In order to understand the details of the anoma- 
lous atomic transport we study the motion of atoms fol- 
lowing their trajectories during the deposition events. 




A. Atomic trajectories in the upper layers 

In order to understand more deeply the mechanism of 
atomic injection of Pt, we plot the atomic positions of a 
crossectional slab cut in the middle of the simulation cell 
(with a slab thickness of 15 A) for the top layer atoms 
during vapor deposition of Pt with the average tq ~ 2.8 
A of the elemental values (upper Fig 2) at K. The out- 
of-plane movement of Al atoms is rather strong during 
injection, while no such movements (or with only much 
smaller amplitude) can be seen when injection is sup- 
pressed by setting rg = 3.1 A (the lower panel of Fig 2). 
We conclude that the presence of the impurity atom is 
necessary for disordering the surface layer locally. 

The positions of the energetic atoms are also shown at 
'--^ K simulation in Fig 2 We can see a chaotic dynamics 
and in certain cases turbulent or circulating atomic mo- 
tions around the equilibrium lattice sites. In certain cases 
the circulating atomic trajectories could also be charac- 
terized by vortices (see the enlarged image of the pattern 
in the right inset of the upper panel of Fig 3). The 
turbulent movement of atoms leads to large lateral and 
out-of-plane amplitudes of atomic vibrations. 

We see no such state of surface atoms in the pure Al 
simulation cell nor in the cases when tq is set to above 
a critical value {vc > 3.1 A). The presence of a Pt atom 
is necessary for the emergence of the surface disordered 
matter. Moreover, the onset of the disordered atomic 
motions seem to be a necessary condition of injection of 
Pt atoms (this will be shown in subsection C). The injec- 
tion of Pt requires the collective motion of few vibrating 
surface atoms. The assistance of at least two transient 
substrate adatoms with large vertical amplitudes seems 
to be a necessary condition for the injective mechanism. 

In the upper panel of the left inset Fig 2 the top view 
of the atomic trajectories of a central vibrating atom 
is shown surrounded by 6 other atoms which suggests 
that the injection of Pt requires the collective motion of 
few vibrating surface Al atoms with a hexagonal sym- 
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FIG. 2. The crossectional view of typical trajectories of 
surface Al atoms induced by the deposition of a Pt atom at 
K (the xz cut of the simulation cell, the scale is in A in the 
axes). The positions of the Al atoms are collected up to 1 ps. 
The first neighbor distance tq — 2.8 A. Inset on the right: 
The enlargement of a typical pattern of atomic dynamics of 
a surface Al atom as a crossectional view. Inset on the left: 
Top view (seen from the (111) surface) of lateral atomic tra- 
jectories in the central region of the surface in the top layer. 
The displaced Al atoms are arranged in a hexagonal lattice. 
Middle Figure: Chaotic trajectories obtained under the same 
conditions given above except that the first neighbor distance 
ro is set to 2.5 A. Initial atomic positions are also shown 
with larger filled circles. Lower Figure: Atomic trajectories 
obtained for ro = 3.1 A 



7 



metry. Interestingly the lateral motion of these atoms 
is anisotropic: the atoms follow an elongated nearly one- 
dimensional trajectories. Another atoms surrounding the 
7 "active" Al atoms are much less affected by the impu- 
rity atom, and the 2nd neighbors are nearly unperturbed 
lattice atoms on the surface. Hence, we can say, that 
when injection occurs atomic vibrations should be cou- 
pled with each other in order to open up a channel (an 
empty surface site) for Pt. 

Tuning ro we also observe the onset of chaotic atomic 
trajectories. In the middle panel of Fig. 2 it can be 
seen that the atomic motions become extremely chaotic 
when ro < 2.5 A. Further decrease in ro give rise to 
the ejection of a pair of Al atoms to the vacuum. This 
kind of an approach raises the question whether at least 
a partly chaos assisted mechanism could be responsible 
for the anomalous atomic transport. The analysis of the 
Lyapunov exponent (that is the measure of the chaotic 
processes) of the system [40] could give the answer for 
this question. Preliminary results indeed indicates that 
the presence of the surface impurity induces atomistic 
chaos on the surface of Al(lll) [59]. These results will 
be published elsewhere. 

The local disordered surface state persists until nearly 
less than a ps. Also, we find the ejection of an Al atom to- 
wards the vacuum (becomes a sputtered atom) when the 
EA EAM potential is used. Hence, there are similiraties 
and differences between the CR and EAM potentials. Us- 
ing the CR potential we find no sputtered Al atoms and 
usually a single Al atom is ejected to the surface while 
with the EA Al and Johnson's cross EAM potential 2- 
4 atoms or more are ejected to the surface and one of 
them is released to the vacuum. Using either the simple 
average of the elemental potentials as a hcteronuclear 
potential or the crosspotcntial obtained by the Johnson's 
scheme [55] we also find injection with the EA EAM po- 
tential. Hence, the injective superdiffusion of Pt is not 
the specific artificial feature of the CR potential but is 
rather a generic feature of the EAM-modcl. The validity 
of the predicted superdiffusive features should be checked 
by more sophisticated ab initio MD approach. This could 
be an apparent work to be done in the future. Although, 
ab initio calculations can be carried out for a substrate 
with the number of atoms up to ~ 100 with present day 
supercomputers, therefore the obtained results could also 
be treated with care. 



B. Acceleration towards the surface 

The deposition of impurity atoms over the substrate 
surface leads to the local acceleration of the particle and 
to the impact to the surface with few eVs kinetic en- 
ergy even if initiated by zero velocity. Hence the phe- 
nomenon of acceleration of the deposited impurity atom 
precedes the injection which has already been a well- 
known phenomenon [17,20,57,58]. Recently, Lee et al. 




0,5 1 1,5 2 

time (ps) 

FIG. 3. The kinetic energy of the accelerating deposited 
particle as a function of the time (ps) using the Cleri-Rosato 
tight-binding potential (ro = 2.83 A in the crosspotcnti al. 
The particle is initialized 4 A above the free (111) surface 
of the Al with nearly zero velocity. The plotted points have 
been collected during few dozens of events with randomly 
varied impact points within the {x,y) = -1-1; —1 (A) region 
of the surface. Inset: The peak kinetic energy of the acceler- 
ating deposited particle Pt before the impact to the surface 
(eV) as a function of the first neighbor distance (ro) in the 
heteronucelar potential. The error bars denote standard de- 
viations obtained for few simulation events at each of the ro 
points. 



demonstrated that the low-energy deposition of Ni on 
Al(OOl) leads to serious acceleration of the Ni atoms 
striking the Al(OOl) surface with 3 — 4 eV kinetic energy 
[58]. In another study Wang et al. found the acceleration 
of Au particles over the Ag(llO) surface. Moreover they 
also find the intermixing of the Au atoms with the Ag 
substrate in 11 % of the events in the first 6 ps of the 
simulations [17]. In the rest of the events the Au atoms 
remain above the surface. In certain cases acceleration 
of deposited particles leads to impact mixing (injection, 
e.g. in Pt/Cu) and it has been explained by the large 
negative heat of solution of this system [20]. Although 
this conclusion is based simply on the comparison of few 
diffusion couples which have different heat of solution. 
However, we demonstrate in this article that not heat of 
mixing {AH„i, or heat of solution) is the decisive factor 
of anomalous mixing. 

We find that even if an artificially repulsive hcteronu- 
clear potential is used (^ = in Eq. (3)) impact mixing 
does occur for Pt/Al or Pt/Cu couples. Also, local accel- 
eration does occur with a weakly attractive interaction 
when e.g. a noble gas Ar atom is deposited on Al (in this 
case no mixing occurs) with a peak kinetic energy of ~ 2 
eV. 

The acceleration of the deposited particle can be fol- 
lowed in Fig 3 when the initial velocity is zero. The 
kinetic energy of the accelerating particle is given in eV 
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as a function of time using the CR potential (obtained 
for few tens of events) the peak value of ~ 3 eV is found 
when a Pt atoms is deposited with zero initial velocity. 
The variation of the strength of the heteronuclear inter- 
action in the potential {AH„i) does not affect the speed 
of acceleration and the peak values of the kinetic energy, 
e.g., if we set the preexponential ^ = 10. in Eq. (3), 
we set in a huge value of heat of mixing in the Al-Pt 
potential (and a very strong internuclear attraction with 
a very deep potential energy well in the potential), how- 
ever, we find that the acceleration of Pt is not affected by 
the deeper potential well in the crosspotential. Also, if we 
set in smaller values, or even ^ = could be set in (purely 
repulsive potential, only the Born-Mayer term is nonzero 
given in Eq. (4)), no change in acceleration is observed. 
Hence we conclude that the speed of acceleration is not 
affected by the strength of the Al-Pt potential and the 
generally accepted opinion, that the particle acceleration 
above solid surfaces is due to the large heat of solution 
of the corresponding alloy system, is not supported [20] . 
In other words we find that not the thermodynamic bias, 
built in the cross-interaction (via the depth of the po- 
tential well) drives the acceleration of the particle during 
deposition. 

The variation of tq in the crosspotential does influence, 
however, the speed of acceleration. We find a correlation 
between and the peak kinetic energy of the particle 
(see inset upper Fig. 3). Although there is some scat- 
ter in the data, however, the trend is evident. Adjusting 
the size-mismatch the kinetic energy of the particle im- 
pact can be tuned by varying the relative positions of the 
minima of the potential wells (the ASM). At the average 
ro of the elemental values 2.83 A) the abrupt change 
of the kinetic energy can be seen which reports us that 
the system is increasingly sensitive to the variation of rg 
at the physically realistic values of ASM. Decreasing rg 
we find the amplification of surface instability and the 
enhancement of out-of-plane vibrations. The emergence 
of the injective mixing mechanism could be due to the 
ASM induced surface lattice instability (local impurity 
induced heating up the surface). 

No acceleration of a standing particle is observed when 
the atom is initialized from above 6 A the surface. We 
find the same decay distance of the surface long range 
forces when the cutoff distance of the crosspotential is 
increased above this value. 

Depositing Al atom on Al(lll) acceleration also oc- 
curs (self-atomic acceleration). The most of the depo- 
sition events with various host-substrate couples lead to 
accelerative deposition regardless to the chemical nature 
of the interacting couples. 

The amplification of the out-of-plane vibration of Al 
atoms during deposition gives rise to the oscillation of 
the electronic density p of few surface atoms and to a 
long-range density tail above the surface. The long range 
forces (which initialize and speed up the acceleration) 
could be induced by the long range tail of the embedding 
function as a function of the density p in Eq. (2). Due 
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FIG. 4. The mean square of out-of-plane atomic displace- 
ments (MSD, {-R^)z, A ) of few Al atoms in the topmost layer 
in the injection zone (SxSA )asa function of simulation time 
(ps) during deposition and without deposition (surface wav- 
ing) using the Cleri-Rosato potential. The deposited atom is 
initialized 4.6 A above the surface. Inset: The crossectional 
view of the uppermost layer ai t — 0.95 ps with the approach- 
ing Pt atom (together with the transient Al atoms). 



to the out-of-plane instability of Al the electron density 
decays slowly above the surface and which gives rise to 
slowly damping F{p). This could lead to the enhance- 
ment of the impurity-surface interaction. 



C. The enhancement of surface disordering 

In order to point out the correlation between the vi- 
brational amplitude of the surface atoms and the acceler- 
ation of the deposited Pt atom we plot in Fig 4 the mean 
square of atomic out-of-plane displacements (i?^)z of few 
surface Al atoms in the middle of the simulation cell on 
the surface (using a 5 x 5 A area) as a function of the 
time (ps). In Fig 5 we plot the lateral {R^)xy in order 
to demonstrate the enhancement of lateral movements of 
the surface Al atoms in the " active" region when the in- 
jection takes place. This, together with Fig 5 confirms 
that the abrupt activation of a surface disordering mech- 
anism is necessary for the injection of the Pt atom. In Fig 
4 the abrupt increase of {R^)z at ~ 1 ps has been found 
which indicates the appearance of transient Al adatoms 
on the surface (see inset Fig 4). This transient process 
coincides with the injection of Pt and catalyzes intermix- 
ing in a similar way as proposed by Haftel et al. [19]. 

The injection of the impurity particle starts with the 
acceleration of the Pt atom. In turn the approaching 
particle induces the amplification of the out-of-plane vi- 
bration of few Al atoms which further accelerates the 
impurity particle towards the surface. Finally, at a criti- 
cal proximity to the surface {dptAi ^ 3 A) strong lateral 
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FIG. 5. The mean square of lateral atomic displacements 
(MSDa;i/,A ) of few Al atoms in the topmost layer in the injec- 
tion zone (5 X 5 A ) as a function of simulation time (ps) dur- 
ing deposition and without deposition using the Cleri-Rosato 
potential. The deposited atom is initialized 4.6 A above the 
surface. 



components to the surface atomic displacements sets in 
(the appearance of the disordered surface state) together 
with the out-of-plane components and which leads to the 
injection of the particle. Without the surface instabil- 
ity no injection occurs, the particle becomes an adatom. 
The concerted motion of few surface atoms opens up a 
channel for injection. 

The amplification of the out-of-plane vibration of Al 
atoms during deposition gives rise to the oscillation of 
the electronic density p of few surface atoms and to a 
long-range density tail above the surface. The long range 
forces (which initialize and speed up the acceleration) 
arise from the long range tail of the embedding function 
as a function of the density p in Eq. (2). Due to the 
out-of-plane instability of Al the electron density decays 
slowly above the surface and which gives rise to slowly 
damping F{p). This could lead to the enhancement of 
the impurity-surface interaction. 

Concerning the generalization of the superdiffusivc fea- 
tures of intermixing in Pt/Al, we can say that injection 
occurs in many other transition metal/Al systems. How- 
ever, in certain cases (e.g. Cu/Al) we find that the in- 
jection of the deposited particles appears after the de- 
position of dozens of atoms. In these cases, the slight 
decrease of ro leads also to ultrafast injection. We see 
similar features in the case of the Ni/Au system [42]. As 
already mentioned above, injection also akes place in the 
Pt/Cu system in a very similar way as in Pt/Al. It might 
be the case that Pt shows superdiffusivc features even in 
more media (substrate). It is also interesting, what hap- 
pens if we invert the systems, e.g. the deposition of Al 
on Pt leads to the lack of intcrdiffusion during the time 
scale we can reach. However, if we set ro = 2.1 A, the 
deposited Al atom can be forced to interdiffuse to the Pt 



phase [42]. Above this value no injection occurs. It must 
also be noted that in general we find a transition from 
superdiffusion towards normal diffusion, e.g. when the 
first neighbor distance is tuned. Finally we should em- 
phasize that the observed anomalous atomic transport 
processes and phenomena (injection with superdiffusion, 
surface assisted acceleration, oscillatory deposit-surface 
interaction, disordered surface matter) might not be the 
artifact of the employed interaction potentials nor the 
PARCAS code. The agreement with experimental vapor 
sputter deposition results [21,24] as well as the accurate 
parameterization of the employed Al many body poten- 
tials [15,43,52,54] and the straightforward usage of them 
for the anharmonic effects of Al [27,57] provide solid basis 
for sharing these results with the wider community. 

An important question remains to be resolved is 
whether a specific state of matter is found or the size- 
mismatched host-impurity interaction induced surface 
disordering can be considered as a local low-temperature 
premclting phenomenon. Anyhow, the impurity induced 
low-temperature disordering of the surface layers includ- 
ing few dozens of atoms is an unexpected short lived lo- 
cal state of the surface and calls for further verifications. 
Although its lifetime is short (less than a ps), however, 
this localized spatial and temporal surface state occurs 
at each impact events. This state of matter resembles a 
cascade event, however, its appearance does not depend 
on the impact energy of the impurity particle. 



V. CONCLUSIONS 

In this paper we studied the deposition of Pt on 

Al ( 1 1 1 ) . An oscillatory adsorbatc-surface interaction and 
ultrafast injective inter-layer atomic transport of the de- 
posited vapor atoms has been found for the Pt/Al couple 
using two different types of atomic interaction potentials. 
The soft landing of the deposit induces the disordering 
of the local region of the surface Al atoms. The disor- 
dered state persists up to ^ 0.5 ps which covers the tran- 
sient out-of-plane vibration of few surface Al atoms. The 
anomalous inter-layer transport of the impurity atom to 
the substrate is assisted by the transient out-of-plane mo- 
tion of few surface Al atoms. This collective exchange 
mechanism works even at very low temperatures slightly 
above K and calls for experimental verification. 

Contrary to the general belief that chemical and ther- 
modynamic forces govern particle deposition and inter- 
diffusion no effect of thermodynamic bias has been found 
on acceleration towards the substrate's surface and on in- 
termixing of the deposited particle. Instead we point out 
the role of atomic size-mismatch (ASM) in particle accel- 
eration and in superdiffusivc features. In particular we 
find that intermixing can be tuned by a system param- 
eter which can be given as the ratio of the atomic size 
parameter of the cross-interaction term to the substrate's 
lattice constant (the ASM). Moreover, the chaotic nature 
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of atomic trajectories as well as the lattice instability of 
the surface can also be tuned by ASM. 
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